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When is a quantum heat engine quantum? 
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Quantum thermodynamics studies quantum effects in thermal machines. But when is a heat 
engine, which cyclically interacts with external reservoirs that unavoidably destroy its quantum 
coherence, really quantum? We here use the Leggett-Garg inequality to assess the nonclassical 
properties of a single two-level Otto engine. We provide the complete phase diagram characterizing 
the quantumness of the engine as a function of its parameters and identify three distinct phases. 

We further derive an explicit expression for the transition temperature. 
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The study of thermal machines has been a cornerstone 
of thermodynamics since its origins. The performance 
of these devices is usually characterized by two key pa¬ 
rameters: efficiency and power [T]. While attention has 
for a long time solely focused on idealized infinite-time 
processes—with maximum efficiency but zero power— 
research in the last decades has shifted towards more 
realistic hnite-time cycles—with reduced efficiency and 
non-vanishing power mm- Macroscopic internal com¬ 
bustion engines, refrigerators and heat pumps are notable 
examples of classical thermal machines whose finite-time 
dynamics follow the laws of classical physics m- Quan¬ 
tum heat engines, on the other hand, are microscopic mo¬ 
tors that are described by dynamical equations of motion 
that obey the laws of quantum mechanics [1] . They have 
been extensively studied over the last hfty years [5HTB]. 
and both quantum coherence m and quantum correla¬ 
tions [m [19] have been theoretically shown to be able 
to boost their performance. Meanwhile, concrete pro¬ 
posals to experimentally build such quantum machines 
have been put forward using, for instance, trapped ions 
[23 EH and nanomechanical systems [22ll23]. However, 
contrary to most quantum applications that aim at per¬ 
fectly shielding systems from their environments , 

thermal machines, by their very nature, cyclically inter¬ 
act with heat reservoirs that unavoidably destroy their 
quantum coherence. A crucial issue that therefore needs 
to be addressed, for both theoretical and practical rea¬ 
sons, is how to identify the quantumness of a heat engine. 

Characterizing nonclassicality is of fundamental im¬ 
portance in many areas of quantum physics, from quan¬ 
tum computation [IZIEH] to quantum biology [221 ISO] ■ 
Assessing the nonclassical properties of a given state or 
dynamics is a nontrivial task, however, and, more than 
a century after the birth of quantum theory, the bor¬ 
der between classical and quantum worlds remains fuzzy 
[STHSS]. A common approach to identify quantum be¬ 
havior is to impose classical constraints that are violated 
by quantum mechanics [84j . The assumptions of realism 
and locality thus lead to Bell’s inequality [5S|, while those 
of macroscopic realism and non-invasive measurements to 
the Leggett-Garg inequality (LGI) [36j. These two results 
allow to test whether a system has stronger-than-classical 
spatial or temporal correlations, respectively. The advan¬ 


tage of the LGI is that it applies to a single particle and 
not to a pair of particles as required by Bell’s inequality. 
It has therefore become a valuable theoretical and experi¬ 
mental tool to detect quantumness of individual systems 
m- In the last few years, experimental violations of 
the LGI have been observed in an increasing number of 
single systems, including a superconducting qubit [38j . 
a spin in a diamond defect center [39j . a nuclear spin 
snuni, a photon [321 US] ; a phosphorus impurity in sili¬ 
con |44| , and a single diffusing atom [45] . It is worth men¬ 
tioning that the last two experiments have implemented 
ideal negative measurements, as originally discussed by 
Leggett and Garg. In the following, we employ the LGI 
to assess the quantumness of a paradigmatic quantum 
thermal machine: a single two-level Otto engine jEHUj. 
Our investigations reveal the existence of three different 
regimes which we summarize in a phase diagram as a 
function of the parameters of the engine (see Figs. 3 and 
4 below). We further derive an explicit formula for the 
transition temperature. We finally show that trying to 
operate a thermal machine faster, with the idea of keep¬ 
ing it coherent, actually leads to incoherent dynamics, 
owing to constraints imposed by thermodynamics. 

Quantum Otto engine. We consider a quantum Otto 
engine for a single two-level system with time-dependent 
frequency ujt and Hamiltonian H = where is 

the usual Pauli operator (we set h = ks = 1). The quan¬ 
tum Otto motor is a generalization of the familiar four- 
stroke car engine and its thermodynamic cycle consists 
of the following isentropic and isochoric steps, as shown 
in Fig. 1: (1) Heating: the two-level system is weakly 
coupled to a hot reservoir at temperature Th during time 
Tfi, while its frequency is kept hxed, (2) Expansion: the 
system is isolated and its frequency is unitarily changed 
from UJ 2 to oJi during time ri, (3) Cooling: the system 
weakly interacts with a cold reservoir at temperature Tc 
during time Tc, while its frequency is again held constant, 
(4) Compression: the frequency is unitarily brought back 
to its initial value 002 during time T 2 . 

Work and heat along the different steps may be deter¬ 
mined from infinitesimal variations of the average energy 
of the system, E = (H) = ujP, where P = {<Jz)l‘2^ is the 
polarization. We write accordingly dE = Pduj uidP = 
SW SQ and associate work 6W with changes of fre- 
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FIG. 1. (color online) Four branches of the quantum Otto 
cycle for a two-level system: (1) constant frequency heating 
during Th, (2) unitary expansion during ri, (3) constant fre¬ 
quency cooling during Tc, (4) unitary compression during T 2 . 


quency and heat 6Q with changes of polarization. The 
Otto engine exchanges energy in the form of work during 
compression/expansion phases and in the form of heat 
during heating and cooling. In order to compute these 
quantities along the four branches of the cycle, we specify 
the finite-time dynamics of an arbitrary system observ¬ 
able Xf with the Markovian quantum master equation in 
the Heisenberg picture [46l |47] , 


dX 

dt 


.Lu dX 

+ ^+] + [ct_,X] cr+) 

+ y(’^(w,'r) -h 1) (cr+ [X, cr_] -h [<7+,X] a_) , 


( 1 ) 


where a± denote the raising and lowering spin operators 
and n{u},T) = l/[exp(/3a;) — 1] the bosonic thermal oc¬ 
cupation number at inverse temperature /? = 1/T. The 
coupling constant 70 vanishes during the unitary com¬ 
pression/expansion steps 2) and 4) when the engine is not 
interacting with the heat reservoirs; in this limit Eq. Q 
reduces to the standard Heisenberg equation. 

During the finite-time operation of the engine, irre¬ 
versible work is produced along the compression and 
expansion steps owing to nonadiabatic transitions. We 
model the corresponding increase of the polarization with 
the help of the formula, Pt = Po + {a/Ti)'^t, (z = 1 , 2 ) [TUb 
na. The phenomenological coefficient a is often referred 
to as internal friction, as its presence reduces the perfor¬ 
mance of the engine [inHni- It is directly related to the 
irreversible entropy production in the usual manner, as 
shown in [iS]- The quantum Otto cycle is thus charac¬ 
terized by the following eight parameters, ti, 2 , Th^c, 
‘Th,ci 7o and a. However, not all choices of parameter 
values are physically meaningful as two important condi¬ 
tions need to be fulfilled in order to have a proper heat 
engine: the total work produced by the machine should 
be positive and the cycle closed. Thermodynamics there¬ 
fore restricts the allowed range of the above parameters, 


in particular, that of the heating and cooling times Th^c 
and of the internal friction coefficient a. 

The optimal heating and cooling times were deter¬ 
mined in Ref. | 12 | using the method of Lagrange multipli¬ 
ers. Computing work and heat along the four branches 
of the Otto cycle, assuming, for concreteness, a linear 
driving protocol, uit = cot + loq, the total work produced 
by the engine in its steady state is found to be [ 12 ] (the 
derivation is summarized in [48| for completeness). 


W = —(uj2 — Wl) 


^peq 


(I-a;)(l-y) 


(I - xy) 


( 2 ) 


(I - xy) \Ti T2 


- + -)(i-y) 


-b cr'wi 1 — + — 
T2 Ti 


where = P^« - P^^"? = - tanh(w 2 /( 2 r,i ))/2 -b 

tanh(a;i/( 2 Tc ))/2 is the difference between the equilib¬ 
rium polarizations after complete thermalization with 
the hot and cold reservoirs, x = exp(—ycTc) and y = 
exp(- 7 /ir;i) with -^h,c = 7 o coth(w/(2T?i_c))- Optimizing 
the work with the constraint of a fixed total cycle 
time T = Ti -b T 2 -b T/i -b Tc leads to the Lagrange function, 


C{x, y,X) = WP\\T-\ - InxH - Iny — ri— T 2 ) , 

\ Ic Ih ) 

(3) 

with Lagrange multiplier A. Equating the partial deriva¬ 
tives of £(a;, y. A) with respect to x and y to zero and set¬ 
ting for simplicity 7zi = 7c = 7 eventually yields (keeping 
the frequencies u )\^2 and times ti _2 constant) the optimal 
hot and cold thermalization times [I 2 j . 


1 
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( 4 ) 

( 5 ) 


for an engine with vanishing positive work—this is the 
minimal requirement for a thermal machine to be a heat 
engine. In that case, the total cycle duration r is the 
minimum time needed to obtain a nonnegative work out¬ 
put. This minimum time vanishes in the limit a to zero 
[I2j . The two quantities R and Xmax are defined as, 

R - <^^^1 (l/p + V'^ 2 ) _ AP^g - cr^/r2 

~ Auj(AP<^i-Ta'^/ri)’ AP®?-b cr^/n ' ^ 


The constraint of a closed cycle, that is, of a finite cy¬ 
cle duration r, leads to an additional condition on the 
internal friction coefficient a that reads m, 

2 

— < AP®«. (7) 

T2 


Heating and cooling times t/j c diverge when the above 
bound is saturated, indicating that large nonadiabatic 
effects, generated by the finite-time driving, result in in¬ 
creased, and eventually infinite, thermalization times. 

Expressions (4) and (7) form the basis of our investi¬ 
gation of the quantum nature of the Otto engine. Since 
















3 



FIG. 2 . (color online) Leggett-Garg function K2,{t), Eq. ( 10 1, 
a function of time for four different values of the damping 
constant 7 . The shaded area indicates the classically forbid¬ 
den region K3 > 1 . Above the critical value, 7/2 ~ 012, there 
are no violations of the LGI ([8| for any value of t. 


I? a) 




□ LGI violated ■ ■ ■ LGI not violated 


coherence is preserved along the unitary compression and 
expansion steps, our strategy is to look for violations of 
the LGI during the heating phase, where decoherence is 
expected to be the most adverse. Our analysis may be 
easily extended to the cooling phase as well. 

Leggett-Garg inequality. Under the classical assump¬ 
tions of macroscopic realism and noninvasive measurabil¬ 
ity, the two-time correlation functions Cij of a dichotomic 
observable, Q = ±1, measured at three distinct times U, 
i = (1, 2, 3), satisfy the inequality |55] . 

= C21 + C32 — Gsi < 1 . ( 8 ) 

Quantum mechanics violates the above inequality. A 
value of the Leggett-Garg function K 3 above one is there¬ 
fore a clear signature of nonclassical behavior [37]. 

In order to assess the quantumness of the Otto engine, 
we evaluate the symmetrized two-time correlation func¬ 
tions Cij of the observable for the dynamics given by 
the master equatio n ([l| ) with the help of the quantum 
regression theorem |46| 7 We find |i5] . 

Cij = ^{Wx{ti),tJx{tj)}) = cos [oj2{U - tj)]. 

( 9 ) 

By further considering equally spaced measurement 
times with time separation t, as commonly done m, 
we obtain the Leggett-Garg function. 


K^{t) = 2e 2 ‘cos(w 2 t) —e cos ( 2 W 2 Q • (10) 

The function K^{t) is shown in Fig. 2 for different values 
of the damping parameter 7 . We observe that the LGI ([^ 
is violated for unitary time evolution with 7 = 0 . More¬ 
over, the amplitude of the violations decreases with in¬ 
creasing 7 until the critical value, 7/2 ~ a;2, at which the 


FIG. 3 . (color online) Phase diagram characterizing the quan¬ 
tumness of the Otto engine as a function of the hot tempera¬ 
ture Th and the internal friction coefficient a = aI t2. a) The 
engine behaves nonclassically when the LGI (j^ is violated, 
i.e., when the quantum time Tq is larger than the heating time 
Th- b) We identify three regimes: i) single quantum-classical 
transition for small d, ii) multiple transitions in an interme¬ 
diate domain, iii) classical dynamics for large a. Parameters 
are aii — 10, 6l) 2 = 20, n = 0.01, T2 = 0.1 and Tc = 1. 


time scale of the coherent system dynamics coincides with 
that of the reservoir induced decoherence, is reached. Be¬ 
yond that point there is no violation of the LGI for any 
value of t and the incoherent evolution imposed by the 
hot reservoir prevails. In the sequel, we set the first mea¬ 
surement at the beginning of the heating phase and the 
third measurement at the end, that is, we choose 2 t = Th- 
We characterize the transition from coherent to inco¬ 
herent dynamics by introducing a time Tg defined as the 
largest possible time for which the LGI can be violated: 


Tq = ma.x{2t\K3{t) > 1 } . 


( 11 ) 


The factor 2 comes from the fact that the three mea¬ 
surements span a total time interval of 2t. Owing to 
the oscillating feature of the Leggett-Garg function (10 1 , 


the quantum time Tg is a step function; it decreases with 
increasing reservoir temperature, as expected, see Fig. 3. 

Results and discussion. The behavior of the engine 
depends decisively on the relationship between heating 
time Tfi and quantum time Tg, that is, on the one hand, 
on the internal dynamics of the machine, as dictated by 
thermodynamics, and, on the other hand, on the deco¬ 
herence process induced by the coupling to the reservoir. 
The dynamics of the Otto motor is nonclassical only if the 
quantum time is larger than the heating time, Tg > t^, so 
that the LGI may be violated. As seen from Eqs. (4) and 
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FIG. 4 . (color online) Phase diagram characterizing the quan¬ 
tumness of the Otto engine as a function of the hot and cold 
temperatures and Tc- The engine displays an analogous 
behavior as in Fig. 3 , the cold temperature Tc having a similar 
effect as the internal friction coefficient u. Same parameters 
as in Fig. 3 with a — 2 / 3 . 


(6), the heating time depends on the two temperatures 
Th and Tc and on the internal friction coefficient cr; we 
will use in the following the reduced quantity a = uI t 2 - 

Our main results are summarized in the phase diagram 
shown in Fig. 3. For a fixed cold temperature Tc, we iden¬ 
tify three distinct regimes: i) ct < ac\- for tr smaller than 
a critical value trci) the behavior of the engine is quan¬ 
tum for temperatures Th <Tq = a; 2 /( 2 coth“^( 2 a; 2 / 7 o)); 
owing to the almost vertical shape of Tq, the threshold 
temperature Tq is mostly independent of a in that do¬ 
main; ii) (Tci < (T < d'c 2 ' for intermediate values of a, 
the engine appears quantum in some temperature inter¬ 
vals and classical in some others. This corresponds to 
a transition regime where the boundary between classi¬ 
cal and quantum worlds is blurry; hi) 0-02 < d: for a 
larger than a second critical value (tc 2 , the dynamics of 
the Otto engine is classical for all temperatures Th- Note 
that the tendency of Th to diverge for larger values of a 
is clearly visible in the upper left corner of Fig. 3a). 

The phase diagram in Fig. 3 reveals some subtle inter¬ 
play between thermodynamics and quantum physics. We 
first remark that in the low a phase, the threshold tem¬ 
perature Tq is equivalent to the condition 7/2 ~ 0^2 for a 
two-level system coupled to a single reservoir at temper¬ 
ature T 2 . Here the internal dynamics of the engine which 
is set by thermodynamics does not play any role, and the 
quantumness of the Otto cycle is solely controlled by the 
decoherence process generated by the reservoir. The sit¬ 
uation changes dramatically in the two other phases as 
the heating time Th increases significantly with increasing 
a. In the high a regime, thermodynamics thus imposes 
a purely classical behavior, although decoherence alone 
would predict the existence of a quantum sector for low 
enough temperatures. In other words, the classical or 


quantum nature of the engine is not just governed by the 
coupling to the reservoirs. A second, nonintuitive, ob¬ 
servation is that trying to operate the engine faster in 
order to keep it coherent is bound to fail: shortening the 
compression/expansion steps indeed results in increased 
nonadiabatic excitations induced by a finite a [49], and, 
in turn, to longer thermalization times. The latter even¬ 
tually lead to stronger decoherence and, consequently, 
to incoherent dynamics. The rule ’the faster, the bet¬ 
ter’ which might be true for most quantum applications, 
hence does not apply to quantum heat engines. 

Figure 4 finally presents the influence of the cold reser¬ 
voir temperature Tc (for fixed a) on the quantumness of 
the Otto engine. We note that it is largely similar to that 
of the internal friction a displayed in Fig. 3, as both pa¬ 
rameters have an equivalent role in exciting the system. 
For Tc below the critical value Tci , it does not affect the 
behavior of the engine which is controlled by the same 
threshold temperature Tq given above (the discussion in 
the previous section was based in this Tc regime). In an 
intermediate domain, the dynamics of the engine alter¬ 
nates between quantum and classical as Th is increased. 
Ultimately, for Tc larger that the critical value Tc 2 , the 
LGI is never violated and the engine is purely classical. 

Conclusions. We have used the Leggett-Garg inequal¬ 
ity to completely characterize the quantum properties of 
a two-level Otto engine. We have obtained a phase dia¬ 
gram that allows to identify the parameter regimes where 
the engine behaves nonclassically, and where, therefore, 
quantum resources in the form of coherence or correla¬ 
tions, might be successfully exploited. We have further 
shown that trying to run a thermal machine faster, with 
the hope of beating decoherence, is not a winning strat¬ 
egy, as thermodynamics will make the dynamics classical. 
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APPENDIX 


A. Relation to entropy production 


We here establish a connection between the internal 
friction coefficient a and the irreversible entropy produc¬ 
tion E. We begin by evaluating the heat exchanged dur¬ 
ing the two isochoric steps. Since the frequency is con¬ 
stant, we have in any time interval [ 11 ,^ 2 ], 

Q= [\,P,dt = u;{Pt,-Pt,). (12) 

Jti 

As a result, we find for the heating and cooling steps, 


Qh = W2 (Pb — Pa) 


(13) 


(pB-PA + <y^ 


\ 

['^2 


where A, B,C,D denote the four corners of the thermo¬ 
dynamic cycle m- The entropy variation of the engine 
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vanishes for a full cycle. The overall entropy change is 
then given by the entropy difference in the two reservoirs: 


AS=- 



Qh 




{Pb - Pa) + 




(14) 



The difference in polarization Pb — Pa is further m, 


Pb-Pa^ 


(l-a;)(l-j/) cr^(l - y)(^ + ^) 


l-xy 


1 - xy 


( 15 ) 

and we can thus cast the entropy variation in the form, 


AS = 


+ 


( 

L02 

Uc 

Th 

( UJl 

UJ2 

Uc 

Th 

Wicr^ / 

' 1 

Tc 1 

J2 




(l-a;)(l-y) 


l-xy 


l-xy 


(16) 


1 

Tl 


We here identify three sources of irreversibility: partial 
thermalization {x,y vanish after complete equilibration), 
the internal friction a, and the intrinsic irreversibility of 
the engine given by the term {uj 2 /Th — uJi(Tc). The ef¬ 
ficiency of the Otto engine is indeed only equal to the 
Carnot efficiency, when uJilui 2 = Tc/T^. In that qua¬ 


sistatic limit, the first two terms in Eq. (16) vanish and 
the total entropy change simplify to. 


AS = 


ujia 


1 

T 2 


1 

n 


ojia 

TcT 2 


0 J 2 cr 

ThTl 


(17) 


We recognize the usual expression for the entropy pro¬ 
duction in the long-time limit, AS = 'SijTi |50] with 
Sc = ui\a^ jTc and = uj 2 <j'^/T h- A finite friction coef¬ 
ficient cr therefore leads to a finite entropy production. 


B. Total work produced by the engine 


For completeness, we summarize in this section the 
derivation of the work produced by the engine [H]. The 
adiabatic work during either compression/expansion is. 




/ Pujdt = (ujf — uJi) 
Jo 



(18) 


where LOij denotes the initial and final frequencies uJi 2 - 
Additionally, the irreversible work associated with the 
nonadiabatic driving (and the corresponding increase in 
polarization) along these branches is given by, 





ujPdt = 


iOf + 

2 Ti ' 


(19) 


The total work done during a compression or expansion 
step is therefore the sum, 


W‘°' = (u, - oj,) 




ujf +uji a 


( 20 ) 


= (ujf - UJ^) Po H- ujf. 

'Pi 


The work produced by the engine during one cycle, 
Eq. (2), is the sum over the compression and expansion 
steps (with proper time matching) [T^ . 


C. Correlation functions of the two-level system 


The two-time correlation functions (9) of the two-level 
system may be obtained in the following way. We intro¬ 
duce the vector A having for components the three Pauli 
operators and the unit operator. The latter quantities 
form a complete set of observables for the two-level sys¬ 
tem. The master equation (1) can then be rewritten as a 
matrix differential equation, dA/dt = M ■ A, oi explicitly. 


/Ox\ 


(-1 

—CO 

0 

(Jy 

- 

(jJ 

1 

2 

0 

<Xz 

~ 

0 

0 

-7 

\l) 


\ 0 

0 

0 


0 

0 

-7o 

0 


\ 

((Jx\ 




Cfz 

/ 

\l) 


( 21 ) 


According to the quantum regression theorem [45], the 
equation of motion for the two-time correlation functions 
is the same as that for the operators, and we thus have, 

^{d{t)A{t + r)) = M{d{t)A{t + r)). (22) 

dr 

Since we are only interested in the correlation functions 
of the operator we can restrict ourselves to the upper 
left 2x2 submatrix of M. We hnd. 


dr 


C = 


-i 

^ -i 


C(r), 


with the correlation vector 
C(r) = 

The solution to Eq. 


(XxipCTxit + t)) 
a^{t)ay{t + t)) 


\ ( e. 2 cos(a;T) —e 2 sin(cLiT) 

\ e 2 sin(a;Tj e 2 cos(a;Tj 


is given by, 

2 l 
2 

jt 


with the initial condition, 

{(Xx(Xx){t) 


C[t = 0) = 


^(JxO'y 


){t) 


(23) 

(24) 

C(r = 0), 

(25) 

(26) 


m ■ 


The last equality follows from the algebraic properties of 
the Pauli matrices. The symmetrized correlation func¬ 
tion is equal to the real part of the above correlation 
function and reads, 

1 T 

Cij =-{[(Jx{t + T),ax{t)]) = e 2'^cos(a;r). (27) 
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